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With a purpose of constructing a robust evolution system against numerical instability for in- 
tegrating the Einstein equations, we propose a new formulation by adjusting the ADM evolution 
equations with constraints. We apply an adjusting method proposed by Fiske (2004) which uses 
the norm of the constraints, C^ . One of the advantages of this method is that the effective signa- 
ture of adjusted terms (Lagrange multipliers) for constraint-damping evolution is pre-determined. 
We demonstrate this fact by showing the eigenvalues of constraint propagation equations. We also 
perform numerical tests of this adjusted evolution system using polarized Gowdy-wave propagation, 
which show robust evolutions against the violation of the constraints than that of the standard 
ADM formulation. 

PACS numbers; 04.25.D- 



I. INTRODUCTION 

The standard way for integrating the Einstein equa- 
tions is to split spacetime into space and time. The 
Arnowitt-Dcscr-Misner (ADM) formulation [l|, Q pro- 
vides the fundamental spacetime decompositions. How- 
ever, it is known that the set of the ADM evolution equa- 
tions is not appropriate for numerical simulations such 
as the coalescences of the binary neutron-stars and/or 
black-holes, which are the main targets of gravitational 
wave sources, and which requires quite long-term time 
integration. 

In order to perform an accurate and stable long-term 
numerical simulation in strong gravitational field, we 
need to modify the ADM evolution equations. This is 
called as the "formulation problem in numerical relativ- 
ity" 

The origin of the formulation problem is the violation 
of constraints, which triggers the blow-up of simulations. 
The discretization of equations arises truncation errors 
inevitably, so that we have to adjust the evolution sys- 
tem which is robust for error-growing modes. Several 
formulations are suggested and applied; among them, the 
Baurngarte-Shapiro-Shibata-Nakamura (BSSN) formula- 
tion [y, 0], the generalized-harmonic (GH) formulation 
[1, i] and the Kidder-Scheel-Teukolsky (KST) formula- 
tion [lo| are applied widely for the inspiral black-hole bi- 
nary mergers. (Many numerical simulations are reported, 
but we here cite the works [III [3 for applications of the 
BSSN formulation, [H for the GH formulation, and 111 
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for the KST formulation). There are also many other 
formulations which are waiting to be tested |15l - [l9| . 

The current succeeded large-scale numerical simula- 
tions are applying such modern reformulations, but also 
using the "constraint-damping" technique, which is ob- 
tained by adding the constraint terms to evolution equa- 
tions. The additional constraint-damping terms are re- 
ported to be the key implementation in BSSN and GH 
system (e.g. [13, HI). We [l6, 22, 23,] systematically in- 
vestigated how the additional constraint terms change 
the original evolution systems, under the name "ad- 
justed systems" . As we will review in SecEIl monitor- 
ing the stability of the evolution is equivalent to check 
the constraint propagation equations (dynamical equa- 
tions of constraints). Therefore, we proposed to analyze 
the eigenvalues of the constraint propagation equations, 
which can predict the violation of constraints before we 
try actual simulations. 

Based on the same motivation with this "adjusted sys- 
tem" , Fiske [13] proposed an adjustment which uses the 
norm of constraints, C^, which we call the "C^-adjusted 
formulation". He applied this method to the Maxwell 
equations, and reported that this method reduces the 
constraint violations for a certain range of the coefficient. 
An advantage of this C^-adjusted formulation is that the 
effective signature of the coefficients is pre-determined. 
In this article, we apply the C^-adjusted formulation to 
the ADM evolution equations, since the ADM formu- 
lation is the one of the most basic evolution systems in 
general relativity. We show the eigenvalue analysis of the 
constraint propagation of this set, and also demonstrate 
numerical evolutions. 

Before the numerical relativity groups faced the for- 
mulation problem, Detweiler |25| suggested another ad- 
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justment based on the ADM evolution equations. He 
proposed a particular combination of adjustments which 
make the norm of constraints damp down. The story is 
quite similar to this work. However, Detweiler's method 
is restricted with the maximal slicing condition, K = 0, 
and also the behavior except the flat-space is unknown. 
We also show numerical demonstrations of Detweiler's 
evolution equation for a comparison. 

We compare the violations of the constraints between 
the standard ADM, Detweiler's ADM and C^-adjusted 
ADM formulations. We use the polarized Gowdy-wave 
evolution which is one of the comparison test problems 
as is known to the Apples- with- Apples testbeds [26| . 
The models precisely fixed up to the gauge conditions, 
boundary conditions, and technical parameters, therefore 
testbeds are often used for comparison between formula- 
tions H^-ii]. 

The plan of this article is as follows. We review the 
idea of adjusted systems and C^-adjusted formulation 
in Secini We also describe a recipe for analyzing the 
constraint propagation with its eigenvalue analysis which 
we call the constraint amplification factors (CAFs). In 
Sec lHIl we apply the C^-adjusted formulation to the 
ADM equations and show its CAFs. We also review De- 
tweiler's formulation in SecHH We show our numerical 
evolutions in Scc lIVl and we summarize this article in 
Sec|Vl In this article, we only consider the vacuum space- 
time, but the inclusion of matter is straightforward. 

II. THE IDEA OF ADJUSTED SYSTEMS AND 
C^-ADJUSTED SYSTEMS 

A. The idea of adjusted systems 

We review the general procedure of rewriting the evo- 
lution equations which we call adjusted systems [l^ [H, 
[221 . [23I I . Suppose we have dynamical variables which 
evolve along with the evolution equations, 

ftu' = /(u\a,w\---), (2.1) 

and suppose also that the system has the (first class) 
constraint equations, 

C"'(u",ajw'',---) w 0. (2.2) 

We propose to study the properties of the evolution equa- 
tion of C" (which we call the constraint propagation), 

ftC'' =g(C",ft;C",---), (2.3) 

for predicting the violation behavior of constraints, C"^, in 
time evolution. The equation (|2.3|) is theoretically weakly 
zero, i.e. 9(C° w 0, since the system is supposed to be the 
first class. However, the free numerical evolution with the 
discrctizcd grids introduces constraint violation at least 
the level of truncation error, which sometimes grows to 
stop the simulations. The set of the ADM formulation 



has such a disastrous feature even in the Schwarzschild 
spacetime, as was shown in [23j . 

Such features of the constraint propagation equations, 
(|2.3|) . will be changed when we modify the original evo- 
lution equations. Suppose we add the constraint terms 
to the right-hand-side of (|2.ip as 

dtu" = f{u\ d,u\ ■■■)+ F{C'',d,C^, ■■■), (2.4) 

where F^C"' , • • • ) w in principle but not exactly zero in 
numerical evolutions, then p.3|) will also be modified as 

a^c^ = g(c^a,c^•••) + G(c^a,c^•••)• (2.5) 

Therefore we arc able to control StC" by an appropriate 
adjustment F(C"',5jC°, • • ■ ) in (pii]) . There exist vari- 
ous combinations of F(C°, ajC", ■ • • ) in ([2^ . and ah the 
alternative formulations are using this technique. There- 
fore, our goal is to find out a better way of adjusting the 
evolution equations which realizes dtC"" < 0. 

B. The idea of C^-adjusted formulations 

Fiske proposed an adjustment of the evolution 
equations in the way of 

dtu^ = f{u\ d,u\ (^-^^ , (2.6) 

where is positive-definite constant coefficient, and 
is the norm of constraints which is defined as 

= y CaCd^x. The term {5C^/5u3) is the functional 

derivative of with We call the set of ()2.6p with 
(|2.2p as "C^-adjusted formulation" . The associated con- 
straint propagation equation becomes 

<,.c-MC".ftc-,..o-/<i»x(f;)««(^;;). 

(2.7) 

If we set k'-' so as the second term in the RHS of 
(|2.7p becomes dominant than the first term, then 9tC^ 
becomes negative, which indicates that constraint vio- 
lations are expected to decay to zero. Fiske presented 
some numerical examples in the Maxwell system, and 
concluded that this method actually reduces the con- 
straint violations. He also reported that the coefficient 
K*^ has a practical upper limit in order not to crash sim- 
ulations. 



C. The idea of CAFs 

There are many efforts of re-formulation of the Einstein 
equations which make the evolution equations in an ex- 
plicit first-order hyperbolic form (e.g. fiol. [TtI [sol I3H). 
This is motivated by the expectations that the symmetric 
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hyperbolic system has well-posed properties in its Cauchy 
treatment in many systems and that the boundary treat- 
ment can be improved if we know the characteristic speed 
of the system. The advantage of the standard ADM sys- 
tem [H (compared with the original ADM system is 
reported by Frittelli [s^l from the point of the hyperbolic- 
ity of the constraint propagation equations. However, the 
classification of hypcrbolicity(weakly, strongly or sym- 
metric hyperbolic) only uses the characteristic part of 
evolution equations and ignore the rest. Several numeri- 
cal experiments [s^ reported that such a classification 
is not enough to predict the stability of the evolution 
system, especially for highly non-linear system like the 
Einstein equations. 

In order to investigate the stability structure of p.Sp . 
the authors [2^ proposed the constraint amplification 
factors(CAFs). The CAFs are the eigenvalues of the 
coefficient matrix, M""}, (below), which is the Fourier- 
transformed components of the constraint propagation 
equations, dtC"". That is, 

wheve C''{x,t) ^ J C{k,tY' exp{ik ■ x)d^k. (2.8) 

CAFs include all the contributions of the terms, and en- 
able us to check the eigenvalues. If CAFs have negative 
real-part, the constraints are forced to be diminished. 
Therefore, we expect more stable evolution than a sys- 
tem which has CAFs with positive real-part. If CAFs 
have non-zero imaginary-part, the constraints are sup- 
posed to propagate away. Therefore, we expect more 
stable evolution than a system which has CAFs with zero 
imaginary-part. The discussion and examples are shown 
in [SjliH, where several adjusted- ADM systems Q and 
adjusted-BSSN systems [ifll are proposed. 



equations are 

n = + K^- K,jK'^ « 0, (3.3) 
= DjK', - D,K « 0, (3.4) 

where ^^^R is the three-scalar curvature, ^'^^R = 

and K is the trace-part of the extrinsic curvature, K = 

The constraint propagation equations of the Hamilto- 
nian constraint, p.3p . and the momentum constraints, 
(|3.4p . can be written as 

dtV. = P'D.n - 2aDiM' + 2aK'H - A{D,a)M\ 

(3.5) 

dtM^ = -{l/2)aD,n 4- |3'^DlM^ - {D,a)n 

+ (D,l3''')Me + aKM^, (3.6) 
respectively. 

Now we apply C^-adjustment to the ADM formulation, 
which can be written as 

/ (5C^ \ 

daij = (EHI - Kyijmn T , (3.7) 

V 07mn / 

dtK,, = ^ - KK^Jrnn ijj^ j , (3.8) 

where is the norm of the constraints, which we set 

= y" {n^ + Y^M^MJ)d^x, (3.9) 

and both coefficients of K^ijmn , i^Kijmn are supposed 
to be positive definite. We write [SC^/S^mn) and 
[SC"^ /5Kmn) explicitly as (PTj) and (|X2|) in Appendix 

B. Constraint Propagation with C^-adjusted ADM 
formulation 



III. APPLICATION TO THE ADM 
FORMULATION 

A. The standard ADM formulation and 
C^-adjusted ADM formulation 

We start by presenting the standard ADM formula- 
tion of the Einstein equations. The standard ADM 
evolution equations are written as 



KK, 



Djp^, (3.1) 
J, DiDja 
P'DiK.j, (3.2) 



2IUeK' 



where (7^^ , Kij ) are the induced three-metric and the ex- 
trinsic curvature, (a,/?') are the lapse function and the 
shift vector, Di is the covariant derivative associated with 



lij 



and ^^^Rii is the three Ricci tensor. The constraint 



In this subsection, we discuss the constraint propa- 
gation of the C^-adjusted ADM formulation, by giving 
the CAFs on flat background metric. We show CAFs are 
negative real numbers or complex numbers with negative 
real-part. 

The constraint propagation equations, (|3.5|1 and p.6p . 
arc changed due to C^-adjusted terms. The full expres- 
sions of the constraint propagation equations are shown 
as (|B1|) and (|BTT|) in Appendix [B 

If we fix the background is fiat spacetimc, (a — 1, — 
0,7.y = dij,Kij = 0), then CAFs are easily derived. 
For simplicity, we also set Kyij„m = KKijmn = KS^mSjn, 
where k is positive. The Fourier-transformed equations 
of the constraint propagation equations are 



dt 



n 

-(l/2)iA:, 



-2ifcj 



U 



(3.10) 
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The eigenvalues, A, of the coefficient matrix of p.lOp 
are given by solving 

{X + K\k\'^)'^{X'^ + AX + B) = 0, 

where A = 4K\k\^{\k\^ + 1) and B = \k\^ + 16K^\k\'^. 
Therefore, the four eigenvalues are 

(_^|A:|2^_«|fc|2,A+,A_), (3.11) 

where 

A± = -2K\kW\k\^ + 1) ± |fc|V-l+4K2|fc|2(|fc|2-l)2. 

(3.12) 

From the relation of the coefficients with solutions, 

A++A_ = -A<0, and A+A_ = B > 0, (3.13) 

wc find both the real parts of A+ and A_ arc negative. 
Therefore, we see all four eigenvalues are complex num- 
bers with negative real-part or negative real numbers. 

On the other hand, the CAFs of the standard ADM 
formulation on flat background [k = in (|3.1ip ] are re- 
duced to 

(0,0,±i|fc|), (3.14) 

where the real-part of all of the CAFs are zero. There- 
fore the introduction of the C^-adjusted terms to the 
evolution equations changes the constraint propagation 
equations to a self-decay system. 

More precisely, CAFs depend to |fcp if k ^ 0. This 
indicates that adjusted terms affect to reduce high fre- 
quency error-growing modes. Since wc intend not to 
change the original evolution equations drastically by 
adding adjusted terms, we consider only small k. This 
limits the robustness of the system to the low frequency 
error-growing modes. Therefore the system may stop due 
to the low frequency modes, but the longer evolutions are 
expected to be obtained. 

C. Detweiler's ADM formulation 

Wc review Detweiler's ADM formulation [2^ for a com- 
parison with the C^-adjusted ADM formulation and the 
standard ADM formulation. Detweiler proposed an evo- 
lution system in order to ensure the decay of the norm of 
constraints, dtC^ < 0. His system can be treated as one 
of the adjusted ADM systems and the set of evolution 
equations can be written as 

dt%j = ^+LD-,,j, (3.15) 
(9tAy = dSJ) + LDKrj, (3.16) 
where D^ij = —a^jijH, (3.17) 

(3.18) 



where L is a constant. He found that with this par- 
ticular combination of adjustments, the evolution of the 
norm constraints, C^, can be negative definite when we 
apply the maximal slicing condition, K = 0, for fixing 
the lapse function, a. Note that the effectiveness with 
other gauge conditions is remain unknown. The numer- 
ical demonstrations with Detweiler's ADM formulation 
are presented in [H, [l^ , and there we can see the drastic 
improvements for stability. 

The CAFs of Detweiler's ADM formulation on fiat 
background metric are derived as [22| , 

(^-(L/2)|fcp,-(i/2)|fcp, 

- (4L/3)|fc|2 ± ^ |fc|2{-i + (4/9)L2|fc|2}^ , (3.19) 

which indicates the constraints will damp down if i > 0, 
apparently better feature than the standard ADM for- 
mulation. 

IV. NUMERICAL EXAMPLES 

We demonstrate the damping of constraint violations 
in numerical evolutions using the polarized Gowdy-wave 
spacetime, which is one of the standard tests for compar- 
isons of formulations in numerical relativity as is known 
as the Apples- with- Apples tcstbedsjl^. The tests have 
been used by several grou ps and were reported in the 
same manner (e.g. [27h29| |). 

The testbeds provide three tests of the solutions of 
the Einstein equations: gauge-wave, linear-wave, and 
Gowdy-wave tests. Among these tests, we report only 
on the Gowdy-wave test. This is because the other two 
are based on the fiat backgrounds and the violations of 
constraints are already small, so that the differences of 
evolutions between the ADM, C^-adjusted ADM, and 
Detweiler- ADM are indistinguishable. 

A. Gowdy-wave Testbed 

The metric of the polarized Gowdy-wave is given by 

ds^ = t-^/^c^/\-dt^ + dx") + tie^'dy^ + e'^'dz^), 

(4.1) 

where P and A are functions of x and t. The time coordi- 
nate t is chosen such that time increases as the universe 
expands, this metric is singular at t = which corre- 
sponds to the cosmological singularity. 

For simple forms of the solutions, P and A, are given 
by 

P = Jo(27rt) cos(27rx), (4.2) 
A = -27rtJo(27ri) Ji(27rt) cos^(27ra;) 4- 2Tr^t'^[jQ{2nt) 
+ J2(27rt)] - (1/2){(2^)2[J2(2^) + J2(27r)] 
-27rJo(27r)Ji(2^)}, (4.3) 
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where J„ is the Bessel function. 

Following [2^, the new time coordinate r, which satis- 
fies the harmonic slicing, is obtained by coordinate trans- 
formation as 

t{T) = fce^^, (4.4) 

where k and c are arbitrary constants. We also follow 
[26j for choosing these constants fc, c and initial time 
as 

k - 9.67076981276405, c - 0.002119511921460, 

(4.5) 

to = 9.87532058290982 (4.6) 

in such a way that the lapse function in new time coor- 
dinate is unity and t = t at initial time. 

We also use following parameters specified in [2^ , 

• Simulation domain: x G [—0.5, 0.5], y ~ z = 0. 

• Grid: Xn = -0.5 + (n - (l/2))da;, 71 = 1, ■ • • , 100, 
where dx ~ 1/100. 

• Time step: dt ~ 0.25dx. 

• Boundary conditions: Periodic boundary condition 
in x-dircction and planar symmetry in y- and z- 
dircctions. 

• Gauge conditions: the harmonic slicing and = 0. 

• Scheme: second order iterative Crank-Nicholson. 

Our code passed convergence tests with the second-order 
accuracy. 

B. Constraint violations and the damping of the 
violations 

Figure [T] shows the L2 norm of the Hamiltonian con- 
straint and momentum constraints with a function of 
backward time (— t) in the case of the standard ADM 
formulation, (|3.ip - p.2[) . We see the violations of the mo- 
mentum constraints are larger than that of the Hamilto- 
nian constraint at the initial stage, and both grow larger 
with time. The behavior is well-known, and the starting 
point of the formulation problem. 

We, then, compare the evolutions with three formu- 
lations: (a) the standard ADM formulation p.ip - (|3.2p . 
(b) Detweiler's formulation ((XT5]) - ((XTC1) . and (c) the C^- 
adjusted ADM formulation ([X71) - ((X5| . We tuned the 
parameters L in (a), and Hjijmn and HKijmn in (c) 
within the expected ranges from the eigenvalue analy- 
ses. In the formulation (c), we set n^ijmn = t^'y^im^jn 
and KKijmn = KKSimSjn for simplicity, and optimized 
K.y and kk in their positive ranges. We use L = — 10+^'^ 
and {k^, Kk) = (—10^^ °, — 10~^'^) for the plots, since the 
violation of constraints are minimized at t = —1000 for 
those evolutions. Note that the signatures of {k^,kk) 
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FIG. 1. The L2 norm of the Hamiltonian and momentum 
constraints of the Gowdy-wave evolution using the standard 
ADM formulation. We see that the violation of the momen- 
tum constraints is larger initially, and both violations are 
growing with time. 
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FIG. 2. The L2 norm of the constraints, C^, of the polarized 
Gowdy-wave tests with ADM and two types of adjusted for- 
mulations. The vertical axis is the logarithm of the and 
the horizontal axis is backward time. The solid line (a) is of 
the standard ADM formulation. The dot-dashed line (b) is 
the evolution with Detweiler's ADM with L = -10+^'^. The 
dotted line (c) is the C^-adjusted ADM with = -10"*" 
and Kk = — 10~^'*. We see the lines (a) and (c) almost over- 
lap until t — —500, then the case (c) keeps the L2 norm at 
the level < 10"'', while the lines of (a) and (b) monotonically 
grow larger with oscillations. We confirmed this behavior up 
to t ~ -1700. 



and L are reversed from the expected one in Sec. ITlland 
Sec. mil respectively, since we integrate time backward. 

We plot the L2 norms of of these three formula- 
tions in Figure [H We see the constraint violations of 
(a) (the standard ADM formulation) and (b) (Detweiler's 
formulation) grow larger with time, while that of (c)(C^- 
adjusted ADM formulation) almost coincide with (a) un- 
til t = —500, then the violation of (c) begins smaller than 
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FIG. 3. The magnitudes of the adjusted terms in each equa- 
tions for the evolutions shown in Figure [21 The vertical axis 
is the logarithm of the adjusted terms. The horizontal axis 
is backward time. The lines (bl) and (b2) are the adjusted 
terms (|3.17|) and (|3.18|l respectively. The lines (cl) and (c2) 
are the adjusted terms (|A1|I and (|A2[I respectively. We see 
the adjustments in Detweiler-ADM [the lines (bl) and (b2)] 
decrease with time, which indicates that these contributions 
become less effective. 



(a). The L2 violation level of (c), then, keeps its magni- 
tude at most O(10~^), while those of (a) and (b) mono- 
tonically grow larger with oscillations. Figure [5] shows 
up to i = —1000, but we confirmed this behavior up to 
t = -1700. 

Figure [2] tells us that the effects of Detweiler's adjust- 
ment appear at the initial stage, while C^-adjustment 
contributes at the later stage. The time difference can 
be seen also from the magnitudes of adjustment terms in 
each evolution equations, which we show in Figure[3l The 
lines (bl), (b2), (cl), and (c2) are the norms of D^ij in 
(IXT71) . DK^J in (|XTa . SC^/S-i.j in JM]), and 5C'^/5Kij 
in (jA2p . respectively. 

We see that the L2 norms of the adjusted terms of 
Detweiler's ADM formulation, D^ij and Dkij, decrease, 
while that of the C^-adjusted ADM formulation increase. 
If the magnitudes of the adjusted terms are smaller, the 
effects of the constraint damping become small. There- 
fore, the L2 norm of of Detweiler's ADM formulation 
are not damped down in the later stage in Figure [51 

One possible explanation for the weak effect of De- 
tweiler's adjustment in the later stage is the existence 
of the lapse function, a (and a^, a'^), in the adjusted 
terms in p.l7p - p.l8p . The Gowdy-wave testbed is the 
evolution to the initial singularity of the space-time, and 
the lapse function becomes smaller with evolution. Note 
that in previous works 0, [l^] , we see that the constraint 
violations are damped down in the simulation with De- 
tweiler's ADM formulation, where the lapse function, a, 
is adopted by the geodesic condition. 

In Figure m we plotted the magnitude of the original 
terms and the adjusted terms of C^-adjusted ADM for- 
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FIG. 4. Comparison of the magnitude of the original terms 
and the adjusted terms of the C^-adjusted ADM formulation, 
(|3T7l) - (|a!8ll . The lines (c3) and (c4) are the L2 norm of the 
original terms [the evolution equations of Qij and Kij, (|3.ip 
and (|3.2|) ]. respectively. The lines (c5) and (c6) are the L2 
norm of the adjusted terms, which is the second terms of the 
right-hand side of (|3.7|l and (|3.8|l . respectively. We see the 
adjusted terms are "tiny", compared with the original terms. 



mulation; the first and second terms in (|3.7p and p.Sp . 
We find that there is O(102)~O(10^) of differences be- 
tween them. Therefore, we conclude that the adjust- 
ments do not disturb the original ADM formulation, but 
control the violation of the constraints. We may under- 
stand that higher derivative terms in (|Aip and (|A2p work 
as artificial viscosity terms in numerics. 



C. Parameter dependence of the C^-adjusted ADM 
formulation 

There are two parameters, and kk, in the C^- 
adjusted ADM formulation and we next study the sensi- 
tivity of these two on the damping effect to the constraint 
violation. 

Figure [51 shows the dependences on and kk- In 
Figure [5l (A) . we fix kk = and change . In Figure [5l 
(B), we fix = and change kk- In FigurelHl (A) , we see 
that all the simulations stop soon after the damping effect 
appears. On the other hand, in Figure[5l (B), we see that 
the simulations continue with constraint-damping effects. 
These results suggest kk 7^ or k-^ = is essential to 
keep the constraint-damping effects. 

We think the trigger for stopping evolutions in the 
cases of Figure [51 (A) (when kk = 0) is the term 
M.5'^^'^'^{dadbdcdd'H) which appears in the constraint prop- 
agation equation of the Hamiltonian constraint, (|Bip . 
We evaluated and checked each terms and found that 
jj^abcd exponentially grows in time and dominates the 
other terms in (|Bip before the simulation stops. Since 
Hs"'"^'^ is consists of 7*^7™" [see jM]) and jHH], the 
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FIG. 5. Parameter dependence of the C^-adjusted ADM formulation. The vertical axis is the logarithm of the and the 
horizontal axis is backward time. The left panel (A) is the evolutions with kk = and — — 10~^ °, — lO""^", — 10~* °, — 10~^". 
The right panel (B) is the cases with k-^ = Q and kk = —10^^'^, —10^^ '', — 10~^'^, — lO"**'^. In (A), we see that the simulations 
stop soon after the constraint dumping effect appears. In (B), we see that the simulations continue with constraint-damping 
effects. 



time backward integration of Gowdy spacetime makes 
this term disastrous. So that, in this Gowdy testbed, the 
cases ~ Q reduce this trouble and keep the evolution 
with constraint-damping effects. 

The sudden stops of evolutions in Figure [S] (A) can 
be interpreted due to a non-linear growth of "constraint 
shocks" , since the adjusted terms are highly non-linear 
Q. The robustness against a constraint-shock is hard to 
be proved, but the continuous evolution cases in Figure 
[S] (B) may show that a remedial example is available by 
tuning parameters. 

V. SUMMARY 

In order to construct a robust and stable formulation, 
we proposed a new set of evolution equations, which 
we call the C^-adjusted ADM formulation. We applied 
the adjusting method suggested by Fiske [l^j to the 
ADM formulation. We obtained the evolution equations 
as (|Aip - (|A2[) and the constraint propagation equations, 
(|Bip and (|B11[) . and also discussed the constraint prop- 
agation of this system . We analyzed the constraint am- 
plification factors (CAFs) on the flat background, and 
confirmed that all of the CAFs have negative real-part 
which indicate the damping of the constraint violations. 
We, then, performed numerical tests with the polarized 
Gowdy-wave and showed the damping of the constraint 
violations as expected. 

There are two advantages of the C^-adjusted system. 
One is that we can uniquely determine the form of the 
adjustments. The other is that we can specify the effec- 
tive signature of the coefficiencies (Lagrange multipliers) 
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independent on the background. (The term "effective" 
means that the system has the property of the dam pin g 
constraint violations). In the previous our study |23 |. 
we systematically examined several combinations of ad- 
justments to the ADM evolution equations, and discuss 
the effective signature of those Lagrange multipliers using 
CAFs as the guiding principle. However, the C^-adjusted 
idea, (|2.6p . automatically includes this guiding principle. 
We confirm this fact using CAF-analysis on the flat back- 
groimd. 

The C^-adjusted idea is the one of the useful ideas to 
decide the adjustments with theoretical logic. We are 
now applying this idea also to the BSSN formulation 
which will be presented elsewhere near future. 

We performed the simulation with the C^-adjusted 
ADM formulation on the Gowdy-wave spacetime and 
confirmed the effect of the constraint dumping. We inves- 
tigated the parameter dependencies and found that the 
constraint-damping effect does not continue due to one of 
the adjusted terms. We also found that the Detweiler's 
adjustment [1^ is not so effective against constraint vio- 
lations on this spacetime. Up to this moment, we do not 
yet know how to choose the ranges of parameters which 
arc suitable to damp the constraint violations unless the 
simulations are actually performed. 

It would be helpful if there are methods to monitor 
the order of constraint violations and to maintain them 
by tuning the Lagrange multipliers automatically. Such 
an implementation would make numerical relativity more 
friendly to the beginners. Applications of the controlling 
theories in this direction are in progress. 
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Appendix A: The additional C'^-adjusted terms 

The adjusted terms SC^ /5^„m and bC^ j bK„in in (|3.7p - 
(13.81) are written as 



S-fm 

+ 2AfH™"7W' - 2{^eM2^"'''^)M' - 2M2^"''''^{^eM') - M"'M", (Al) 



sc 



2 



2i/4""'H + 2M3^""^M' - 2(9^M4^""^)A^' - 2M4^""^\^eM'), (A2) 



where 



jjmn 2fV^^ -\- i"^) ^'^'^ (3)pme6(3)pr 



eb 



- 2KK""' + 2K'"j K"^ , (A3) 



fjimn _ _^fm(3)pTi _ ^fn(3)pTn _|_ ^mn(3)p« _j_ (3)p' 



71771 £ 



+ 



(3)pmn£ (3)p^nm 



(A4) 
(A5) 
(A6) 



ajn^bn 



(A7) 



+ (A8) 

(A9) 



1777-71 ^/"IJv^ 

4 i 



jL^mn J^i7n7n j^^j^^ j^7nn ^^mn . ^ ]\/f-^^^ " 

^^jmn^ are the same with the appendix of [2^ if (m, n) = 
(n, m). 



Appendix B: The constraint propagation equations of C^-adjusted ADM formulation 

The propagation equation of the Hamiltonian constraint with C^-adjustcd ADM formulation can be written as 

+ MsJ^idbd^M'') + mQa'''''idbdcddM''), (Bl) 



9 



where 



14"'" 



- 2H^^^'\dkd,H2'n + 2i73'""^'=(9fec>rf9ei?3'''')} - 2(afe9,K^„„y)|i73'""''Hl' - H:i^''^'{d,H2'n 

+ i73""'='(9daci?3'^''^)| - 2KK™„.y i/4'""i/4^^ (B2) 

+ H2"''"'{dddcH3'^'^'') + 2i72""'^(a<?aci?3'^°') + Hi^'"'\diHi'^) + H3""''"'{dkHi'^) - Hi^'"'HdtdcH2'"') 

+ 2i73"""(afca,5e7?3'^"')| - 2(9,K^™„y)|-772""'i?2'^'' + 2i72'""'(9ci?3'^"'=) + 2i73""'"i?i'^' 

- 2iJ3'""^"(aci?2*^'') - 2i?3™"^'=(afci72*^'°) + 2i73""^"(ad9ci?3'^''') + 4i^3™"^'=(9fc9ci^3'^'''')| 

- 2(afc9,K^„„„,)|-i73""''i?2^^" + 2i/3'"""(5cif3'^"^)|, (B3) 

ab o J zir mn tt ijab tt mna tt ijb , o tj rnna f ct tt ijbc\ . tj mni f r\ tt iiab\ , rr mnab tt ij 

^ -ZhjmnijS -til J^3 — J^2 ^2 + ^^2 (CcJ^S j + -"2 [OeJ^3 ) + J^3 J^l 

_ H:i"''^''\dcH2'''') - Hj,"''"'\diH2'''') - i73""''"(9fci?2'^'') + i^3"""''(9d9ci^3'^''') 

+ 2i73""''^(9^9ei?3'^'^') + 2i/3'""'^"(9fc9,i/3*^''=) + i73'""'^^(afe9^i?3'^"'')| 

- 2(9fK^„„ij)|i72""^i?3'^'''' - 2i?3""'^"H2'^'' + 4i?3""'^"(9,i73^J'"=) + 2i/^""^'=(9fciJ3'J"^)| 

- 2(afe5,K^™„,,)if3'""''i?3'^°', (B4) 

I TJ- mna TT ijbc tt mnab tj iic 1 ou mnab ( tt ijce\ 1 rj mnai / ^ tt ijbc\ 1 rj mnka { c\ tt iibc\ I 
■-imni]Sn2 113 — J^3 J^2 + ^tlj, [Oe-tls ) + -"3 [Oe-tl3 ) + -"3 {C'kJ^3 )> 

- 4(9fcK^™„y )773""'^''i?3'^'^ (B5) 

Oi.- IT mnab TT ijcd (Tl(l\ 

~Zhjmnij-ti3 J^3 1 (,£>Dj 

= -2a(3)r''6, - 4a^, - k^„„„, |2Hi'""Mi,^^' - 2Hi^^{ddM2a'"') - + 2i?2""'(a£Mia'J') 

- 2i72""'(a,a<iAf2a^^'') + 2i/3'"""'(9fc9£Mi,^^) - 2H:r''''\dkd^ddM2a'"')^ 

- (9fK^™„,j)|2i/2""'Mla'^' - 2H2'^''\ddM2a'"') ~ ff2""'A^^'<^^\ + 4^3'""^'= (^feMi,*^' ) 

- 4i/3'™"=(5fcadM2a'^'^)| - (9,9,K^„„,j) {2i/3""''Afia*^ - 2i/3'""''(9rfM2a'^'') - i/3""''A1^'^^\} 

- KKmm, ^2H^''' M:,a'' - 2^4™" (9f Af4a^^0 | , (B7) 



= -2k^ 
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- H2"'"''M^S\ ~ H2"''^''M'S\ + 2H3''""''\diMiJ^) + 2H:i""'^\dkMia'^) - 2Hj,""'^\diddM2a'^'^) 



_ 2H3^-''''M^6\ - 2i/3""''A^'<5^a| + 2{dkd,K^„,mj)W'""''M2a'"' + 2KK,nmjHr"M4a'^'', (B8) 

_ 2i/3'™'=^(aa/2a^^^) - H3"'^''-M'S\ - Hs'^^'-M'Pa^ + 4(9fe«^„„, )i^3"""'M2a^^^ (B9) 
Hg,"^'^ = 2K^„,™,i?3""''^M2a^^''. (BIO) 

The propagation equation of the momentum constraint with C^-adjusted ADM formulation can be written as 

(Bll) 



where 



- M2j^''\ded,H2'n + A/2a""'(a,a<ia,i/3'^''^)| - 2{deK^^m,)^M2,r''H^'' - A/2a""'(aei?2^^^) 
+ Af2a""'(9d9ci?3'^''^)| - 2KA-™™,|Af3a'""i/4'^ + Af4a""' (9£H4'^' ) | " 2(a,/«K„„y)Af4a'™'i?4'^ 

(B12) 

M2a' = -ia^a" - 2/«^„„„,|-A//ia'""i/2^^'' + 2Mia"'" {dcH^'^'"') + Aha'^'^'Hi'' " A'/2a'""'(9cif2'^'^) 

- Af2a'""'(9fif2'^') + Af2a""'(5daei?3'^''^) + 2A/2a""'(9£a,i73''^'^) | 

- 2(5,K^„„,)|-Af2a""'i?2'^' + 2Af2a™"'(a,i73^^''^)| - 2.^^™™, Af4a""''i?4*-'' , (B13) 
Ml3a''^ = -2AC^„„,y |Afia""i73'^"^ - Af2a""''i?2'^'^ + 2Af2a'""'(c'di?3'^'^'') + Af2a""' (a,i73^^'"=) | 

~ 2(9,/^^„„,,)M2a'""^i^3'^■'^ (B14) 

4[4a'^' = -2AC^„„,,M2a""''i/3^^'^', (B15) 

M5ab = 7mb/3'".a + P'lab,e + aK^ab " K^™™, |2Afi,'"" Afib'^' - 2A/i<,™"(9dA//2h^^'^) - A//i,""M(^5')fc 

+ 2Af2a™"'(9fA^ib'^") - 2A//2a""'(5,9dAf26'^'')| - (9fAC^™„y)|2Af2a'""'Afifc^^ - 2Af2a'""'(a<iAf26'^'') 

- A//2a""'X(^ <5')fc| - 2«;K™™, |A/3a""Af3b '^' ^ A/3a"" (9f Af4h '^'' ) + Aha"'"' (deAhb'' ) 

- M4a"'"\deddAUb'''')^ - 2(9,K^^„„,j)|Af4a""^Af36'^' - A/4a""'(9dA/4h*^'^)|, (B16) 



11 

- 2/«Kmn,j|-M3a""A//46^^^ + A//4a""^M3fa'^' - Af4a™"^(9dM46^^'*) - Af4a""' (a,Af46^^'^) | 

+ 2(a,«Kmn,J)A^4a""'M46^^^ (B17) 
My.b^'* = 2K^„„„M2,""^Af26^^'' + 2K^f„„,Af4a'""^M4b^^''*. (B18) 
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